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Simulating a quantum system is more efficient on a quantum computer than on a classical com- 
puter. The time required for solving the Schrodinger equation to obtain molecular energies has been 
demonstrated to scale polynomially with system size on a quantum computer, in contrast to the well- 
known result of exponential scaling on a classical computer. In this paper, we present a quantum 
algorithm to obtain the energy spectrum of molecular systems based on the multi-configurational 
self-consistent field (MCSCF) wave function. By using a MCSCF wave function as the initial guess, 
the excited states are accessible; Entire potential energy surfaces of molecules can be studied more 
efficiently than if the simpler Hartree-Fock guess was employed. We show that a small increase 
of the MCSCF space can dramatically increase the success probability of the quantum algorithm, 
even in regions of the potential energy surface that are far from the equilibrium geometry. For 
the treatment of larger systems, a mult i- reference configuration interaction approach is suggested. 
We demonstrate that such an algorithm can be used to obtain the energy spectrum of the water 
molecule. 

PACS numbers: 



I. INTRODUCTION 

Since the discovery of a polynomial quantum algorithm 
for factorization [l[ , other quantum algorithms that pro- 
vide exponential speedup over their classical counter- 
parts have been found. Examples in diverse areas in- 
clude the computation of approximations to the Jones 
polynomial Q and certain instances of the hidden sub- 
group problem Q. Feynman observed that simulating 
a quantum system might be more efficient on a quan- 
tum computer than on a classical computer 4]. Fur- 
ther work by others has born out this early sugges- 
tion @, 1, 0, 1, i, M, EH- Although a quantum com- 
puter to carry out the calculations that we propose is not 
currently experimentally realizable, many recent develop- 
ments in quantum information technology [l2l Il3l . [Til Hq 
continue to get closer towards the implementation of such 
a device. 

In quantum chemistry, where molecular quantum sys- 
tems are simulated on a classical computer, one is re- 
stricted to employ a finite basis to span the formally in- 
finite Hilbert space that would describe the electronic 
structure of a molecular system. The full configuration 
interaction (FCI) method diagonalizes the molecular 
Hamiltonian to provide solutions to the electronic struc- 
ture problem that are exact within this basis. FCI scales 
exponentially with respect to the size of the molecular 
system studied and therefore is restricted to the treat- 
ment of small diatomic and triatomic systems (l~7j . Re- 
cently, a quantum algorithm for the solution of the FCI 
problem in polynomial time was proposed by Aspuru- 
Guzik et ai. [9J. This algorithm employed the HF wave 



function as a reference for further treatment of the cor- 
relation effects by the FCI Hamiltonian on the quantum 
computer. The excited states of molecular systems are 
difficult to resolve by employing the HF wave function 
as an initial trial state. The main reason for this dif- 
ficulty is due to the fact that contributions from sev- 
eral configuration state functions (CSF) must be consid- 
ered if one is seeking a reasonable overlap of the trial 
state with the exact wave function. In the quantum 
chemical study of molecular systems, people are often 
interested in computing molecular properties, such as 
the energy of the ground state and a few low-lying ex- 
cited states, i.e., in study of the spectroscopic properties 
of molecules. In such cases, an FCI calculation might 
become too expensive even for a quantum computer 
for some large systems. A multi-reference configura- 
tion interaction (MRCI)-truncated Cl-calculation based 
on an a multiconfigurational self-consistent field (MC- 
SCF) wave function can sometimes provide results within 
chemical accuracy, but with much less computational 
work than FCI due to the smaller Hilbert space asso- 
ciated with the calculation. It is difficult to describe 
various regions of molecular potential energy surfaces, 
sometimes even qualitatively correct, by using a single 
reference determinant. Many reference determinants or 
configuration state functions are often required for the 
description of bond-dissociation regions. 

In this paper, we suggest a quantum algorithm to ob- 
tain energy eigenvalues of a MRCI wave function of a 
molecular system using the MCSCF wave function as 
initial input to a quantum computer. We show that by 
improving the quality of the trial wave function, the pro- 
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posed algorithm yields substantially higher success prob- 
abilities than by employing the HF wave function. The 
use of a MCSCF wave function simultaneously reduces 
the amount of quantum computing resources needed and 
extends the range of reliable quantum computations to 
excited states and treacherous regions of the potential en- 
ergy surface. Simulating a chemical system with a quan- 
tum computer requires the mapping of the Fock space 
of the MCSCF wave function to the Hilbert space of the 
quantum bits (qubits) of a quantum computer. We intro- 
duce a more compact mapping technique for molecules 
by employing symmetry properties. This approach re- 
duces the computational resources for representing the 
wave function on a quantum computer and avoids the 
state-crossing problem. 

The structure of this work is as follows. In Sec. |TT] 
we will review the implementation of the FCI scheme on 
a quantum computer. Sec. IIIII describes the properties 
of the MCSCF wave function. In Sec. IIVI we describe 
a quantum algorithm for using MCSCF trial wave func- 
tions in a FCI quantum algorithm. In Sec. [V] we discuss 
numerical evidence for the feasibility of this scheme as 
applied to calculations for the water molecule. We final- 
ize with a conclusions section. 



II. IMPLEMENTATION OF CI SCHEME ON A 
QUANTUM COMPUTER 

A closed quantum system in the non-relativistic limit 
can be described by its Schrodinger equation (atomic 
units are used), 



E = E n . 



M fr. 



(i) 



Feit [18| and coworkers suggested a method to solve the 
Schrodinger equation based on the spectral properties of 
the solutions to the time-dependent Schrodinger equa- 
tion. Its solution can be expressed as a linear superposi- 
tions of eigenfunctions of the Hamiltonian, 



tp(r, t) = A n u n (r)exp(-iE n t) 



(2) 



where the function u n (r) satisfies the equation Hu n — 
E n u n . The method requires a numerical solution of 
\4>(r,t)) and the correlation function P(t): 

P(t) = (V(r, 0)| V(r,i)) = J r(r,0)^(r,t)dr, (3) 

where |^(r, 0)) is the wave function at t = 0. P(t) can 
then be expressed as 



P(t) =^K| 2 expH£n*), 



(4) 



P(E) = J2\A n \ 2 5(E~E n ). 



(5) 



A scheme similar to the one proposed by Feit can 
be implemented on a quantum computer. Abrams and 
Lloyd suggested finding eigenvalues and eigenvectors 
using a quantum phase estimation technique. Eigen- 
functions of the Hamiltonian are also eigenfunctions of 
the unitary time-evolution operator, U(t) = exp(— iHt), 
whose eigenvalues can be expressed as a phase factor. 
A quantum Fourier transform [12] (QFT) is used to re- 
trieve the phase in a binary expansion and thus obtain 
the eigenenergy. This scheme has been proposed to sim- 
ulate quantum systems, especially Fermion systems, on 
a quantum computer 0, [2Q] . If the Hamiltonian 

can be decomposed by means of a split-operator tech- 
nique @, H, [ll|, [l8| , the quantum computational cost is 
polynomial, it can provide an exponential speed increase 
over its classical counterpart. 

The details of the algorithm proceeds as follows 0, [H| : 
First, one must prepare two quantum registers, one is the 
index register composed of m qubits, which are used as 
control qubits and to perform a QFT operation. Another 
register of n qubits is the target register that is used to 
represent the wave function of the system. The index 
register is initially prepared in the zero state |0). The 
quantum bits of the index register are entangled with 
successive binary powers of the unitary evolution opera- 
tor on the target register. After the time-evolution of the 
target register, the index register encodes an eigenvalue 
of the time evolution operator U of the target system as 
a phase represented in a binary notation. By performing 
a QFT, the phase, and therefore the eigenvalue of the 
system can be obtained. 

The algorithm begins by initializing the quantum com- 
puter into the state: 



l*o> = \0M- 



(6) 



Performing a tt/2 rotation on each qubit in the index 
register results on the state 



Af-l 



(7) 



3=0 



where M = 2"\ By performing a series of controlled-C/ 
operations on this state, it is transformed into: 



M-l 



*2> = i £ u j \m). 



(8) 



3=0 



The approximate vector \ip) can be written as a sum of 
eigenvectors of U, 



which can be Fourier transformed to display the energy 
spectrum of the system as a set of sharp local maxima at 



(9) 
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where k sums over the dimensionality of the target reg- function is: 
ister. The eigenvalue associated with \<fi k ) is e^ fc , which 
can be written as e 27rilJ fc/ M ; where ujk £ [0,M). Using 
this fact, the state can be rewritten as: 



MCSCF — / , 
K 



D K <f>K, 



(11) 



M-l 

|* 2 > = J>|0 fc )-= J2 e 2 ^ M \j). (10) 

k v j=0 

A QFT performed on the index qubits will reveal the 
phases u>k and thereby the eigenvalues. It requires ~ m 2 
operations [22l |. A polynomial number of trials are re- 
quired to obtain any eigenvalue for which the correspond- 
ing eigenvector is not exponentially small in the initial 
guess. If the initial guess is close to the desired state, 
then only a few trials may be necessary. Once a mea- 
surement is made and an eigenvalue is determined, the 
target register qubits will collapse into the state of the 
corresponding eigenvector. 

Aspuru-Guzik et al. extended the algorithm to the 
study of molecular systems and simplified the algorithm 
by introducing a recursive phase-estimation technique 
that saves the qubits for performing the phase estima- 
tion. They also introduced an adiabatic state prepa- 
ration (ASP) technique for obtaining molecular ground 
states. They demonstrated that such algorithms can be 
applied to problems of chemical interest using modest 
numbers of quantum bits. 



III. MCSCF WAVE FUNCTION 

In the previous quantum computing for quantum 
chemistry work [9|, a HF wave function was used as the 
initial trial wave function. The HF method [23| repre- 
sents the wave function as a single Slater determinant. 
In most cases, the HF wave function by itself is not suf- 
ficiently accurate to generate useful chemical predictions 
such as relative energies of products and reactants, and 
therefore a correlated calculation is necessary. Often, the 
HF wave function is not a good initial guess to the exact 
wave function of the system, especially for the excited 
states calculations, in which contributions from several 
Slater determinants must be considered, even for a qual- 
itatively correct description. If a number of electron con- 
figurations are relatively close in energy (i.e. degenerate 
or near-degenerate), then the HF approximation is par- 
ticularly poor. This is the usual case when one explores 
regions of avoided crossings (or anti-crossings), molecules 
close to the dissociation limit, in the limit of large system 
size, or in the study of a chemical reaction path [24j In 
such cases, it is more appropriate to describe the system 
with more appropriate wave functions in which several 
different electron configurations are taken into account. 

One realization of such a wave function comes from 
MCSCF theory. The general form of an MCSCF wave 



$K = (N\r 1 ^det\l[ ( t> i \, (12) 

i<ZK 

which is a linear combination of several electron config- 
uration state functions (CSF). Each CSF differs in how 
the electrons are distributed between the molecular or- 
bitals (MOs), c/) i . For a particular system, the CSFs can 
be chosen based on physical consideration of the system. 
The MOs are usually expanded in a basis of atomic or- 
bitals (AOs), x„. To obtain a MCSCF wave function, 
both the configuration expansion coefficients Dk and 
the MO expansion coefficients C^i are variationally opti- 
mized. Hence the optimized vector is the best approxima- 
tion to the exact wave function of the system in a specific 
parameter space. For a given set of orbital and config- 
uration parameters, even in a small variational space, 
the MCSCF wave function can give a much better ap- 
proximation than the HF wave function. A truncated CI 
based on an MCSCF wave function, the so-called multi- 
reference CI method, normally, give better results than a 
CI using a HF wave function as a reference, when small 
Hilbert spaces are involved. The trade-off between the 
MRCI approach and the FCI approach is that chemical 
intuition is involved in selecting the appropriate CSFs for 
constructing the CI expansion. 

As mentioned above, computational resource require- 
ments are significantly less for any reasonable MCSCF 
calculation than for an FCI calculation in the same or- 
bital space. Simple combinatorial arguments show that 
/ 2M \ 

there are I „ I possible Slater determinants formed 

from M molecular orbitals and N electrons, and al- 
though their number can be reduced by space- and spin- 
symmetry considerations, the growth in the number of 
determinants with system size remains exponential. Even 
a well-constructed algorithm that uses an iterative pro- 
cess for a subset of the roots (such as that by Lanzcos 
or Davidson) [25| will have CPU requirements that scale 
roughly as the square of the number of determinants. 
Moreover, the storage requirements scale as the number 
of determinants. Consequently, the FCI problem scales 
exponentially with system size. 

In contrast, MCSCF uses an iterative process to obtain 
an optimal (in the variational sense) space of specific size, 
e.g., 8 electrons in 12 orbitals. Even if one adopts the 
most costly (but simply definable) MCSCF calculation, 
the so-called Complete Active Space SCF (CASSCF), 
and the resulting number of determinants is given by the 
same formula as for full CI, the CASSCF space is a tiny 
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fraction of the full CI space. The number of iterations 
required to determine the optimal space is usually on 
the order of 10. The non-CI part of an MCSCF calcula- 
tion is typically dominated by the transformation of elec- 
tron repulsion integrals from the atomic basis in which 
they are calculated to the molecular orbital basis, which 
scales as M 5 . In fact, in real MCSCF calculations, the 
integral transformation step is often the limiting step. 
It is worth noting that more complex MCSCF calcula- 
tions can be made possible by the use of the macrocon- 
figuation approach [2^ |. which can reduce the number 
of determinants in an MCSCF to a polynomial number 
even for larger orbital spaces. This essentially guarantees 
that a physically meaningful and mathematically robust 
MCSCF calculation will be i ntegral bound, and therefore 
scale as M 5 with system size [23, [28 1 for systems of tens of 



atoms. Asymptotically, MCSCF has an exponential cost 
as well and therefore a quantum computer still provides 
an exponential speedup for this method. 

In the MCSCF method, several states can be cal- 
culated simultaneously through a state-averaged ap- 
proach [HHil. For the n-th MCSCF CI root, the energy 
function can be written as, 



(^n\H\A, 



(14) 



A more general energy-like function can be constructed 
by use of weighting vector [3l[ , 



E 



(15) 



where Wi is the weight for state i. So, if we are interested 
in a few evenly or non-evenly weighted states, the MO 
expansion coefficients are optimized for all these states. 
By diagonalizing the one-particle density matrix, we can 
obtain the occupation numbers in Fock space for each 
state. This will be used as initial guess and map to the 
qubits on a quantum computer in the quantum algorithm 
proposed in this work. 



qubits. The direct mapping has the advantage of yield- 
ing a simple Trotter expansion in terms of a polynomial 
number of second-quantized Fermion operators. 

The compact mapping considers the restriction of the 
multiplicity of the system and reduces the number of 
qubits to represent the wave function to a Hilbert space 
where all the quantum states of the target register cor- 
respond to valid electronic configurations within a given 
spin symmetry. The challenge of employing the compact 
mapping to general quantum systems is that the repre- 
sentation of the time-evolution operator may involve a 
larger number of non-local quantum gates. 

Here we introduce a more compact mapping tech- 
nique, which considers the symmetry restriction of the 
molecules. The electronic states can be categorized into 
different irreducible representation of their point group. 
The subspace associated with a particular irreducible 
representation can be mapped to the Hilbert space of 
the target register. This results in considerable savings 
in the number of qubits required to represent the wave 
function. Since there is no interaction between states 
that belong to different irreducible representations, this 
mapping technique can aid in solving certain cases of the 
state crossing problem [l6[ . 

For the proposed scheme, the wave function of the de- 
sired state is implemented as the initial input to the phase 
estimation algorithm using the MCSCF approach. The 
approximation to the exact wave function of the i-th state 
I*,-), is \^f ICSCF }. The probability of observing the ex- 
act i-th state is |(*i|*f CSCF )\ 2 . Since the MCSCF wave 
function provides a much better approximation to the 
ground state wave function of the system than does the 
HF wave function, and also provides a better description 
of excited states than a Koopmans' theorem estimate [l6| 
from a HF wave function, the probability of obtaining the 
correct energy of the system in the phase estimation pro- 
cedure is higher for MCSCF wave functions than for HF 
wave functions. 

The first step for the quantum algorithm involves the 
preparation of a MCSCF calculation for N states that are 
of interest for the system. The MCSCF wave function for 
the state of interest are used as the initial guess for the 
trial wave function: 



IV. IMPLEMENTATION OF A GENERAL CI 
ALGORITHM BASED ON AN MCSCF WAVE 
FUNCTION ON A QUANTUM COMPUTER 

The first step for the proposed simulation algorithm is 
to map the wave function of the system to the state of the 
target register. In quantum chemistry basis set methods, 
many-particle molecular wave functions are represented 
in terms of a single-particle basis expanded in terms of 
atomic orbitals and a many-particle basis expanded in 
terms of Slater determinants or CSFs. In the direct map- 
ping [9|], each qubit represents the fermionic occupation 
state of a particular atomic orbital. The Fock space of the 
molecular system is mapped to the Hilbert space of the 



I*) 



(16) 



where is the MCSCF wave function for the n-th 

state. The next step is to map the MCSCF wave function 
for the n-th state as the initial input to the quantum com- 
puter. This will be prepared using a state-preparation 
algorithm. General state preparation is a hard problem, 
but generally the MCSCF wave function contains a poly- 
nomial number of non-zero terms in the Hilbert space, 
and therefore may be prepared efficiently [13]. Feeding 
the MCSCF wave function into the phase estimation al- 
gorithm as initial guess, the eigen-energies of the corre- 
sponding CI state can be retrieved. 
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An MCSCF vector can be expanded as follows: 

hO = a?) 
k 

where \ip k ) is the eigenvector of the CI matrix. |c„| 2 = 
K^nl^n)! 2 i s th e probability of obtaining the eigenvector 
|V> n )- A ^1 vector for the n-th state can be written as: 

ivo = w:) + 10 = W> + hO + io, (is) 

where is the part of the CI vector in the model 

space, which is used to construct the MCSCF wave func- 
tion; j^) is the part of the CI vector in the space exter- 
nal to the model space; \ifit ev ) i s the deviation of MCSCF 
wave function |^/) n ) from \ip™}, the projection of the CI 
vector in the model space. Then we have: 

(^K) = l+(^|Vr) + (^K)- (19) 

The vectors in model space and external space are or- 
thogonal, = (CIC) = °- We can see that if 
the deviation vector goes to 0, the overlap of the MC- 
SCF vector with the CI vector is one, the algorithm will 
be deterministic. 



V. APPLICATION TO THE WATER 
MOLECULE 

We have performed a quantum simulation for the 
ground state and the first singlet excited state of the 
water molecule using the cc-pVDZ basis set [33| . For the 
ground state, considering the C2V symmetry of the water 
molecule, the HF wave function of water is: 

(la 1 ) 2 (2a 1 ) 2 (16 2 ) 2 (3a 1 ) 2 (16 1 ) 2 . (20) 

We consider a complete active space (CAS) type MC- 
SCF method: the first two a\ orbitals are frozen, the 
active space consists of 3ai — 6ai orbitals, lb\, and I&2 
and 2£>2 orbitals. The MRCI is performed using the same 
model space but considering the single and double exci- 
tations to the external space. The MCSCF space con- 
tains 152 CSFs. The CI space contains 13872 CSFs, here 
logl 3872 — 13.76, so 14 qubits are required to represent 
the CI wave function on a quantum computer. The ge- 
ometry used in the calculation is near the equilibrium 
geometry (R = 1.8435a and LHOH = 110.57). We 
varied both OH bonds from 0.5 to 10 times of the equi- 
librium distance simultaneously, keeping the sym- 
metry, R = aRo,a = 0.5 — 10. The success probabil- 
ity of the quantum algorithm for using HF and MC- 
SCF wave function as initial input | (Wf F \^f T ) | 2 and 
^mgsgf^ci^ ^ are gnown m p ig 1 By following 

the stretch coordinate, we observe that the success prob- 
ability for using MCSCF wave function as initial guess 
is very high (> 0.9) through the stretching, while the 
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FIG. 1: Success probability (P = | < ty HF \-<H CI > | 2 and 
P=\< ^mosof^oi > |2) of uging HF and MCSCF wave 

function as the initial guess. Black line is for HF wave func- 
tion, red line is for the MCSCF wave function of the ground 
state, green line is the MCSCF wave function for the excited 
state. The system is the water molecule, where a = R/Ro is 
the ratio between the stretched bond length R and the bond 
length near equilibrium distance Ro = 1.8435ao. 
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FIG. 2: Success probability (p>mcscf^ci > ^ of uging a 
few CSFs as the initial guess. Black line is for the ground state 
and red line is for the excited state. The system is the water 
molecule, where a — R/Ro is the ratio between the stretched 
bond length R and the bond length near equilibrium distance 
Ro. 



success probability for using HF wave function as initial 
guess decreases very fast as the OH bond is stretched. 
We can still obtain high probability of success by using 
just a few CSFs instead of all 152 CSFs in the MCSCF 
model space. In Fig. 2, we show the success probability 
for both the ground state and excited states using 6 and 
8 CSFs respectively. With a relatively small number of 
CSFs one can have a reasonable overlap with the desired 
state. 

We further studied the performance of the method for 
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TABLE I: Results for the first singlet excited state of the water molecule using the phase estimation algorithm. The MCSCF 
and HF wave function are used as initial guesses. The FCI energy is -83.464130 a.u., the exact energy for using MCSCF wave 
function is -83.449186 a.u. and for using HF wave function is -83.443206 a.u. 



Digits (qubits) 


Energy (MCSCF) 


Energy (HF) 


2 


-83 ± 2.07 x 10 1 


-83 ± 2.07 x 10 1 


8 


-83.6386 ± 5.04 x 10" 1 


-83.6386 ± 5.04 x 10 _1 


16 


-83.4486 ± 1.28 x 10 -3 


-83.4435 ± 1.27 x 10 -3 


24 


-83.44919786 ± 4.95 x 10~ 6 


-83.44318182 ±4.95 x 10 -6 



excited states. We explored the first excited state of the 
water molecule at the equilibrium geometry using the 
STO-3G basis set [34]. The first two a\ orbitals were 
frozen. The model space for the MCSCF is a complete ac- 
tive space that includes the 3ai, 4ai, l&i and 1&2 orbitals. 
The MRCI calculation uses the same model space, but 
considers the single and double excitation to the external 
space. 

We use the scheme introduced by Parker and Ple- 
nio [HI to implement the QFT. This method is known as 
the measured quantum Fourier Transform (mQFT) ap- 
proach. In this scheme, only one control qubit is used, 
more qubits are saved for representing the wave function. 
The mQFT approach is based on the fact that the gates 
within the Fourier transform are applied sequentially on 
the qubits. Thus instead of performing the entire trans- 
form and then making measurements on all control qubits 
afterwards, one can apply the single qubit operation to 
the first qubit and then measure it. The operations con- 
trolled by this first qubit are then replaced by single qubit 
operations given the result of the measurement on the 
first. The measurement outcome is fed back into the 
quantum calculation and this procedure is recycled till 
all the required binary digits are resolved. The target 
register must remain coherent during the whole proce- 
dure. For more details on this procedure, the reader is 
referred to Ref . [35[ . 

For the excited-state simulation, the CI space is com- 
posed of 18 CSFs, so 5 qubits are required to represent 
the wave function. In Table I, we present the results 
for the calculation of the first excited state of the wa- 
ter molecule in the STO-3G basis using the mQFT algo- 
rithm. The MCSCF wave function and HF wave function 
are used as different initial guesses. The MRCI energies 
are obtained to different digits of accuracy depending on 
the number of ancillary control qubits employed. The er- 
ror bars in the table come from the numbers of the qubits 
in the index register. The more control qubits in the in- 
dex register, the more binary digits can be retrieved. For 
example, if n qubits are used as control qubits, then one 
can only obtain up to n binary digits of accuracy in the 
phase estimation, all the binary digits after these n dig- 
its will be uncertain. Therefore, the error is the same 
regardless of the initial trial state (HF or MCSF) em- 
ployed. The FCI energy is in this case is —83.464130 
a.u. The first singlet excited state energy for the water 
molecule using 24 qubits (E=-83. 44919786 a.u.) is lower, 



even including the error bars, than the exact energy us- 
ing the MCSCF wave function (E=-83.449186 a.u.), this 
is because the error in expansion of the unitary matrix is 
only up to the second order in Trotter expansion Q . 



VI. DISCUSSION AND CONCLUSIONS 

In certain regions of molecular potential energy sur- 
faces, electronic states can cross each other or have low 
gaps, like in the case of avoided crossing regions or at 
the bond-dissociation limit. In these cases, the interac- 
tions between states should be considered simultaneously. 
Consequently, the use of single determinant based meth- 
ods is challenging. 

Using an MCSCF wave function as the initial guess can 
deal with the strong interaction between states straight- 
forwardly. This can avoid possible convergence of the 
state wave function to some undesired and unphysical 
states when the energy gap between these states is small. 

By using the more compact mapping technique, cross- 
ing states that belong to different irreducible representa- 
tions can be addressed separately since there is no inter- 
action between the states. For states in an avoided cross- 
ing region and at the dissociation limit where states are 
near degenerate, since the interaction has been consid- 
ered qualitatively in the MCSCF calculation, the overlap 
of the MCSCF wave functions with the corresponding CI 
wave functions are still large, so that even in such regions 
the probability for the reference states is high. Therefore, 
we conclude that the MCSCF wave function can be used 
as a good initial guess for correlated wave functions us- 
ing quantum computing to explore the whole potential 
energy surfaces for ground and excited states with high 
probability of success. 

Using an HF wave function as the initial guess chooses 

a path H HF — > H, for the evolution from the HF state 

to the CI state. In our scheme, we choose the path 
jjmcscf ^ and the stateg evolve from the MCSCF 

state to the MRCI state. Unlike in the case of HF wave 
function in which the evolution is started from a sin- 
gle element of the CI matrix, the MCSCF wave function 
starts the evolution from a small matrix. This makes the 
evolution safer and faster, especially for a MRCI space. 
From the simulation we can see that by including a few 
CSFs in the initial guess, the success probability can be 
increased from very small to near unity. This idea might 
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be used in developing other quantum algorithms. 

VII. ACKNOWLEDGMENT 

We would like to thank Brian Austin for critical read- 
ing of the paper, Wanyi Jiang, Peter J. Love and Ivan 
Kassal for helpful discussions. We acknowledge the sup- 



port of Purdue Research Foundation and the Harvard 
Faculty of Arts and Sciences for their financial support. 
M. R. H. acknowledges the National Science Founda- 
tion (NSF-0313907) for support of his contribution to 
the research presented herein. A. A. G. thanks the Army 
Research Office (ARO-52944-PH) for their support for 
conducting this research. 



[1] P. Shor, in Proc. 35th Annu. Symp. on the Founda- [18 
Hons of Computer Science (ed. Goldwasser, S.) 124- 
134 (IEEE Computer Society Press, Los Alamitos, [19 
California, 1994). 

[2] D. Aharanov, V. Jones and Z. Landau, Proceedings [20 
of the Annual ACM Symposium on Theory of Com- 
puting, 427-436 (2006). ' [21 

[3] D. Bacon, A. Childs, J±6th Annual IEEE Symposium 

on the Foundations of Computer Science, 469-478 [22 
(2005). 

[4] R. Feynman, Inter. J. Theor. Phys. 21, 467 (1982). 
[5] S. Lloyd, Science 273, 1073 (1996). [23 
[6] C. Zalka, Proc. Roy. Soc. Lond. A 454, 313 (1998). 
[7] D. Abrams, S. Lloyd, Phys. Rev. Lett. 83, 5162 

(1999) . [24 
[8] D. Lidar, H. Wang, Phsy. Rev. E., 59, 2429 (1999) 
[9] A. Aspuru-Guzik, A. D. Dutoi, P. J. Love, M. Head- [25 

Gordon, Science 309, 1704 (2005). [26 
[10] F. Remade, R. Levine, PNAS 101, 12091 (2004). 
[11] I. Kassal, S. P. Jor dan, P. J. Love, M . Mohse ni, and [27 

A. Aspuru-Guzik, \http://arxiv. org/ abs/ 0801 . 2986\ 

Submitted, (2008) [28 
[12] H. Plantenberg, P. C. de Groot, C. J. P. M. Harmans, 

J. E. Mooij, Nature 447, 846 (2007). [29 
[13] L. Aolita, L. Davidovich, K. Kim, H. Haifner, Phys. 

Rev. A 75, 052337 (2007). [30 
[14] J. R. Petta, A. C. Johnson, J. M. Taylor, E. A. Laird, 

A. Yacoby, M. D. Lukin, CM. Marcus, M. P. Han- [31 

son, A. C. Gossard, Science, 309, 2180 (2005). 
[15] B. P. Lanyon, T. J. Weinhold, N. K. Langford, M. [32 

Barbicri, D. F. V. James, A. Gilchrist, and A. G. 

White, Phys. Rev. Lett. 99, 250505 (2007) [33 
[16] T. Helgaker, P. Jorgenson, J. Olsen, Molecular [34 

Electronic- Structure Theory, John Wiley & Sons Ltd 

(2000) . [35] 
[17] L. Thogersen, J. Olsen, Chem. Phys. Lett. 393, 36 

(2004). 



M. Feit, J. Fleck, A. Steiger, J. Comp. Phys. 47, 412 
(1982). 

D. Abrams, S. Lloyd, Phys. Rev. Lett. 79, 2586 
(1997) 



S. Wiesner, http://arxiv.org/abs/quant-ph/9603028 
(1996). 

B. C. Travaglione and G. J. Milburn, Phys. Rev. A. 
63, 032301 (2001). 

M. Nielsen, I. Chuang, Quantum Computation and 
Quantum Communication, Cambridge Univ. Press, 
Cambridge (2000). 

A. Szabo, N. Ostlund, Modern Quantum Chemistry: 
Introduction to advanced Electronic Structure The- 
ory, 108, McGraw-Hill (1989). 

M. W. Schmidt, M. S. Gordon, Annu. Rev. Phys. 
Chem. 49, 233 (1998). 

E. R. Davidson, J. Comp. Phys. 17, 87 (1975). 

Y. Khait, J. Song, M. R. Hoffmann, Int. J. Quantum 

Chem. 99, 210 (2004). 

T. Tsuchiya, B. O. Roos, Mol. Phys. 104, 1123 
(2006). 

G. Barone, A. Silvestri, B. O. Roos, Phys. Chem. 
Chem. Phys. 7, 2126 (2005). 

H. Werner, W. Meyer, J. Chem. Phys. 74, 5802 
(1981). 

R. Differenderfer, D. Yarkony, J. Phys. Chem. 86, 
5098 (1982). 

M. R. Hoffmann, C. Sherrill, M. Leininger, H. Schae- 

fer III, Chem. Phys. Lett. 355, 183 (2002). 

A. Soklakov, R. Schack, Phys. Rev. A. 73, 012307 

(2006). 

T. H. Dunning, Jr., J. Chem. Phys. 90, 1007 (1989). 
W. J. Hehre, R. F. Stewart and J. A. Pople, J. Chem. 
Phys. 51, 2657 (1969). 

S. Parker, M. Plenio, Phys. Rev. Lett. 85, 3049 
(2000). 



